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Abstract 

The gluino contributions to the C ' 7 8 Wilson coefficients for h —> 57 are calculated 
within the unconstrained MSSM. New stringent bounds on the 5 ^ and d ^ 1 mass 
insertion parameters are obtained in the limit in which the SM and SUSY contribu¬ 
tions to 6*7,8 approximately cancel. Such a cancellation can plausibly appear within 
several classes of SUSY breaking models in which the trilinear couplings exhibit a 
factorized structure proportional to the Yukawa matrices. Assuming this cancella¬ 
tion takes place, we perform an analysis of the h —> 57 decay. We show that in a 
supersymmetric world such an alternative is reasonable and it is possible to saturate 
the b —> 57 branching ratio and produce a CP asymmetry of up to 20 %, from only 
the gluino contribution to Cj 8 coefficients. Using photon polarization a LR asym¬ 
metry can be defined that in principle allows for the 67,8 and 67 8 contributions to 
the b —> sj decay to be disentangled. In this scenario no constraints on the “sign of 
H” can be derived. 

1 Introduction 

The precision measurements of the inclusive radiative decay B —> X s ^ provides an im¬ 
portant benchmark for the Standard Model (SM) and New Physics (NP) models at the 
weak-scale, such as low-energy supersymmetric (SUSY) models. In the SM, flavor chang¬ 
ing neutral currents (FCNC) are forbidden at tree level. The first SM contribution to the 
b —> s'y transition appears at one loop level due to the CKM flavor changing structure, 
showing the characteristic Cabibbo suppression. NP contributions to b —> sy typically 
also arise at one loop, and in general can be much larger than the SM contributions if no 
mechanisms for suppressing the new sources of flavor violation exist. 
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Experimentally, the inclusive B —> X s 7 branching ratio has been measured by ALEPH 
|H], BELLE [2] and CLEO 0 resulting in the current experimental weighted average 

BR{B -► X s7 ) exp = (3.23 ± 0.41) x 1(T 4 , (1) 

with new results expected shortly from BABAR and BELLE which could further reduce the 
experimental errors. Squeezing the theoretical uncertainties down to the 10% level has been 
(and still is) a crucial task. The SM theoretical prediction has been the subject of intensive 
theoretical investigation in the past several years. From the original calculation at LO (4J , 
impressive progress in the theoretical precision has been achieved with the completion of 
NLO QCD calculations || 0. || and the addition of several further refinements [[|. [1C| . 
The original complete SM NLO calculation [0 gives the following prediction for ^fz = 
m c /mb = 0.29: 

BR(B -> X s7 ) S M = (3.28 ± 0.33) x 1(T 4 . (2) 

The main source of uncertainty of the previous result is due to NNLO QCD ambiguities. 
In [[IT!] it is shown that using \fz = 0.22 (i.e. the running charm mass instead of the pole 
mass) is more justifiable and causes an enhancement of about 10% of the b —> s 7 branching 
ratio, leading to the current preferred value: 

BR(B -> X s7 ) SM = (3.73 ± 0.30) x 1(T 4 . (3) 


Although these theoretical uncertainties can be addressed only with a complete NNLO 
calculation, the SM value for the branching ratio is in agreement with the experimental 
measurement within the 1 — 2 cr level. 

The general agreement between the SM theoretical prediction and the experimental re¬ 
sults have provided useful guidelines for constraining the parameter space of models with 
NP present at the electroweak scale, such as the 2HDM and the minimal supersymmetric 
standard model (MSSM). In SUSY models superpartners and charged Higgs loops con¬ 
tribute to b —> s 7 , with contributions that typically rival the SM one in size. To get a 
sense of the typical magnitudes of the SUSY contribution to b —> s 7 , it is illustrative to 
consider the (unphysical) limit of unbroken SUSY but broken electroweak gauge symme¬ 
try, which corresponds to the supersymmetric Higgsino mass parameter fi set to zero, and 
the ratio of Higgs vacuum expectation values tan/5 = v u /vd set to 1. In this limit SM and 
SUSY contributions are identical in size and cancel each other [1^|, due to the usual sign 


difference between boson and fermion loops. Of course, this limit is unphysical: not only 
must SUSY be (softly) broken, but /i = 0 and tan/5 = 1 have been ruled out by direct and 
indirect searches at LEP. 

I 11 the realistic case of softly broken SUSY, the contributions to b —> sy depend strongly 
on the parameters of the SSB Lagrangian, as well as the values of /1 and tan/5. In par¬ 
ticular, as the origin and dynamical mechanism of SUSY breaking are unknown, there 
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is no reason a priori to expect that the soft parameters will be flavor-blind (or violate 
flavor in the same way as the SM). Of course, the kaon system has provided strong FCNC 
constraints for the mixing of the first and second generations which severely limit the pos¬ 
sibility of flavor violation in that sector |L3], [b|]. Note however that the constraints for 
third generation mixings are significantly weaker, with b —> sy providing usually the most 
stringent constraints. 

Nevertheless, for calculational ease one of the following simplified MSSM scenarios have 
often been assumed: 


The SUSY partners are very heavy and their contribution decouples, so that only 
the Higgs sector contributes to b —» sy. In this scenario, as well in general 2HDMs, 
NLO calculations have been performed |8], [15], [16|. Due to coherent contributions 
between SM and Higgs sector, a lower bound on the charged Higgs mass can usually 


be derived [|17j] in this class of models. In the large tan/3 region the two-loop SUSY 
correction to the Higgs vertex can produce quite sizeable modifications and should 
be carefully taken into account []T8| . 


The SUSY partners as well the extra Higgs bosons have masses of order the elec- 
troweak scale, but the only source of flavor violation is in the CKM matrix. This 
scenario, known as minimal flavor violation (MFV), is motivated for example within 
minimal supergravity (mSUGRA) models. MFV scenarios have been studied at LO 
in certain limits at NLO 


20 , 21 , 22 


23 


and including large tan/3 enhanced 
two-loop SUSY contributions [|T^, In this scenario, the b —> sy decay receives a 
contribution from the chargino sector as well as from the charged Higgs sector. To 
avoid overproducing b —► sy, the charged Higgs and chargino loops must cancel to 
a good degree. This cancellation can be achieved for a particular “sign of p” in the 
mSUGRA parameter space] 1 ] which flips the sign of the chargino contribution relative 
to the SM and charged Higgs loops, always interfering constructively. Although this 
cancellation can occur and puts important constraints on the mSUGRA parameter 
space, it is important to note that it is not due to any known symmetry but rather 
should be interpreted, in an certain sense, as a fine-tuning. 


• There are new sources of flavor violation in the soft breaking terms. In this case, addi¬ 
tional SUSY loops involving down-type squarks and gluinos or neutralinos (hereafter 
neglected compared with the gluino loops due to the weaker coupling) contribute to 
b —> sy. It is well known that the gluino contribution can dominate the amplitude 
for such nonminimal SUSY models, both due to the a s /a enhancement with respect 
to the other SM and SUSY contributions, and due to the m,g/mb enhancement from 

1 Specifically the relative sign between the parameters p and A t: and so generally different from the 
“sign of p v relevant in the case of the muon g — 2 MSSM contribution | 25 ||. 
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the chirality flip along the gluino line. Thus in this scenario, which is generally 
noted as the unconstrained MSSM (uMSSM), usually only the gluino contribution 
is discussed. It has been shown m n that the 23-LR off-diagonal entry of the 
down-squark mass matrix is severely constrained by b —> sq measurements to be of 
0( 10~ 2 ). Less stringent bounds can be obtained for the other 23 off-diagonal entries. 
No known symmetry assures that these constraints can be automatically satisfied; 
again this fact could be interpreted at the electroweak scale as a fine-tuning. 


A discussion of the b —> sq process in the general unconstrained MSSM is in principle pos¬ 
sible, but it is necessary to deal with two unavoidable problems: (i) a large number of 
free, essentially unconstrained parameters, and (ii) the need to achieve a quite accurate 
cancellation between the sizeable different contributions (SM, Higgs, chargino/neutralino 
and gluino) to the Wilson coefficient C 7 associated with the Q 7 oc m b sj_,a ,ll, bjiF lw operator 
in such a way that the experimental measurement, which approximately saturated solely 
by the SM result, is satisfied. Moreover, in general MSSM models with nonminimal fla¬ 
vor violation the gluino loop can also contribute significantly to the Wilson coefficient C 7 
associated with the chirality-flipped operator, Q 7 oc as has been recently 

emphasized in the literature [^7], ^3|. However, as the SM, Higgs, and chargino contribu¬ 
tions to C 7 are typically suppressed by a factor of 0(m s /m b ), it is not possible in general 
to achieve a cancellation between the different terms in C 7 and thus a stronger fine-tuning 
has to be imposed. 

However, it has been recently shown |2f| that in many classes of SUSY breaking models 
a particular structure of the soft trilinear couplings A of the soft-breaking Lagrangian can 
be derived which can alleviate these constraints. Writing these couplings as A tJ = AijYij 
(in which Y denotes the fermion Yukawa matrices), the matrices A for the up and down 
sector are given respectively by: 


/i( M ) /R 4- A^’ u 

^ij ^ii ' t/i 


/\(d) __ A L 1 \R,d 
ij ii ' "jb 7 


(4) 


As shown in |29|, this factorization holds quite generally in string models, for example 
in Calabi-Yau models in the large T limit or in Type I models ]3(J, as well as in gauge- 
mediated [pi] and anomaly-mediated models p2l [33], [34]| . If eq 


holds, specific relations 
can be derived for the off-diagonal LR entries in squark mass matrix. In particular, the 
leading contribution to the entries of interest for the b —> sq process are given in the SCKM 
basis as: 
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with V^'ji the rotation matrices for the up and down quark sector from the interaction 
to the mass eigenstate^ From eqs.(|^-§), one can realize hrst that the down-sector LR 
off-diagonal entries are naturally suppressed by a factor of 0 (rrib/m t ) compared with the 
up-squark sector ones due to the particular factorization of the soft trilinear couplings 
given in eq. (f|) . Second, in these classes of models both the 23 and 32 entries are of the 
same order and proportional to the largest mass (up or down). Consequently, in these 
classes of models, 0( ICC 2 3 ) off-diagonal entries in the down-squark sector along with 0(1) 
off-diagonal entries in the up-squark can be considered in some sense as a prediction of 
the underlying fundamental theory^. This fact implies comparable chargino and gluino 
contributions to b —> sy, making the possibility of cancellations between the W and the 
different SUSY contributions to the Q 7 operator less unnatural. The constraints on the 
gluino contribution to Q 7 are simultaneously alleviated. This flavor structure holds in 
essentially all attempts to build string-motivated models of the soft-breaking Lagrangian. 

The structure of the paper is as follows. In section |2|, we briefly summarize the the¬ 
oretical framework for the calculation of the b —► 57 branching ratio at LO and NLO. In 
section j|, we derive useful mass insertion (MI) formulas for the gluino contributions to the 
Wilson coefficients C' 7.8 and C^g.We demonstrate explicitly that in the large tan (3 region, 
a good understanding of these expressions is obtained only by retaining terms in the MI 
expansion through the second order. For fi of the same order as the common squark mass 
parameter and large tan /?, new (previously overlooked) off-diagonal terms become relevant 
in the b —> sy process. We then devote our attention in section [4] to the analysis of the 
gluino contribution to C 7 8 in the general uMSSM. I 11 particular, we ask the question of 
whether the contribution to C 7 alone can saturate the b —> sq branching ratio, assuming 
that the SM and SUSY contributions to C 7 cancel each other to an extent that the effects 
of C 7 are subleading. While this scenario may initially appear to be unnatural, we will 
argue that sufficient cancellations in C 7 do not involve significantly more fine tuning than 
the usual cancellation required in MFV scenarios. With this analysis, we thus provide an 
alternative interpretation of b —> sy which is at least as viable as any supersymmetric one. 
This analysis also provides more general mass insertion bounds on S 23 than those obtained 
recently [^HJ, where the SM (and sometimes Higgs and chargino) contributions to C 7 are 
always retained. As we are generally interested in moderate to large values of tan/? , we 

2 In this notation the CKM matrix is Vckm = (V^V ■ 

3 It is important to note however that the off-diagonal entries of A in the SCKM basis contain terms 
proportional to the products of entries of the left-handed and right-handed quark rotation matrices, which 
are largely unconstrained (except for the CKM constraint for the left-handed up and down quark rotation 
matrices which enter (for example) ^ 23 )- The quark rotation matrices are highly model-dependent. While 
the diagonal entries can in general safely to be taken 0 ( 1 ), it is typically assumed that the off-diagonal 
quark rotation matrices are suppressed by powers of the Cabibbo angle in a way that mirrors the CKM 
matrix (see e.g. ^|). Note though that this assumption is not required, particularly for the right-handed 
quark rotation matrices which enter T 32 which are of particular relevance for this paper. 
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are able to put rather stringent bounds on the mass insertion parameter | R . In section 
|4.3| , we study the branching ratio and CP asymmetry as functions of the SUSY parameter 
space within this scenario, assuming complex off-diagonal Mis. Throughout the paper, to 
avoid EDM constraints we set the relevant reparameterization invariant combinations of 
the flavor-independent phases to zero. Finally in section |4.4| we show that if the photon 
polarization will be measured, it is possible to distinguish such a scenario from the usual 
C-J dominated scenario through the definition of a LR asymmetry. 

Since we are interested in analyzing a supersymmetric world where the one-loop SUSY 
effects are of the same order as the SM loops, we assume relatively light superpartner 
masses. Specifically we choose the gluino mass rhg = 350 GeV and the common diagonal 
down-squark mass rhc = 500 GeV, with the lightest down-squark mass in the 250 — 500 
GeV range. All of the other sfermion masses, as well the chargino and neutralino masses, 
do not enter directly in our analysis and (some of them) can be taken to be reasonably 
light as suggested by [ [T5|| . Motivated by the lower limit on the Higgs boson mass 
(which suggests | cos 2 /?| ~ 1 ) and by the muon g — 2 excess, we focus to some extent on 
moderate to large values of tan/?, though our formulas and much of the analysis hold in 
general. 

2 b —> §7 branching ratio at NLO 

For the purpose of presentation, we summarize the theoretical framework for evaluating 


the b —> sq branching ratio at NLO. A complete and detailed discussion can be found for 
example in [7. |8|. 0. The starting point in the calculation of the B meson decay rates is 
the low-energy effective Hamiltonian, at the bottom mass scale //&: 

H-eff 

AQ 

(9) 

The operators relevant to the b - 

sq process 4 are: 


Q 2 

= s L 'y ll c L c L 'y lx bL , 


Q 7 

= 2 m b s L ai lt 'b R F IJl , , 


Qs 

= -^m b s L a^G;„T a b R . 

(10) 

and their L R chirality counterpart: 


Q 2 

= snnpCRCR'fbR , 


4 This of course depends on the basis chosen; we have chosen the one easiest for our discussion. 
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Qt = ■^^m b s R (T lu 'b L F ltv , 

Os = ^rm.s^G^TM ■ (11) 

The Wilson coefficients C^j 8 are initially evaluated at the electroweak or soft SUSY 
breaking scale, which we generically denote as /x 0 , and then evolved down to the bottom 
mass scale /x 6 . The standard 1 ’! RG equations for the U 7 2 , 7 ,s operators from the electroweak 
scale (nw < m t) to the low-energy scale ji b is given by: 


62 (/ifc) 

= \ (v ™ 62 (fj-w) , 

( 12 ) 

CM/iQ 

= r]^C 7 (fi W ) + | - V™) c 8(jiw) + h ih ai , 

6 i =1 

(13) 

^8 (fib) 

8 

= rj^Cs(iiw)+ J2 h iV ai , 

(14) 


i =1 


The C' 2 7 8 


where r) = a s (/iw)/ce s (^b) and hi, hi and a* are constants (see 0 ] for details), 
coefficients obey the same running as their chirality conjugate counterparts. If the NP 
scale is much higher than m t , the running from /isusy to pw with six quarks should also 
be taken into account (see the first paper of |I 8 |). The coefficient C 2 is dominated by a SM 
tree-level diagram and is normalized such that C 2 (l-iw) = 1- Its chirality conjugate, C 2 , has 
no SM contribution at tree level and can thus be safely set to zero. The NP contributions 
to C 2 and C 2 appear at one-loop order and are negligible. The Wilson coefficients C 7 and 
C 7 are the only coefficients that contribute directly to the b —> sy branching ratio at the 
lowest QCD order (cc°). These coefficients receive contributions both from the SM and NP 
at one-loop order. The coefficients Cg and C' s receive one-loop SM and NP contributions 
through the same types of diagrams as C 7 and C 7 , but with the external photon line 
substituted by a gluon line. When the QCD running from the matching scale po to fi b is 
performed, these different coefficients mix, as shown in eqs. (|T2l - |T4D , so that the “effective” 
low-energy coefficients 62 , 7,8 (^&) receives contributions from different operators. 

The b —> 57 branching ratio is usually defined by normalizing it to the semileptonic 
b —> c e~ v e branching ratio, giving: 


BR(B 


X s ^)\e i >{1-S)E^ 


= BR(B -> X c ev) 


6 a 

nf{z) 


v tb v* 


V r 


cb 


K(S, z) 


(15) 


5 In a recent paper |27| it has been pointed out that the gluino contribution (and the same argument 
holds also for the chargino and neutralino contributions) is the sum of two different pieces, one proportional 
to the bottom mass and one proportional to the gluino mass, which have a different RG evolution. We 
have found that at LO, this is equivalent to the usual SM evolution once the running bottom mass m b (no) 
is used instead of the pole mass in the Ciffio) WC. 
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Here f(z) is a phase space function and should be calculated for on-shell masses, namely 
y/z = m c /m b = 0.29. 5 is the experimental photon detection threshold, which for com¬ 
parison between experimental data and theoretical prediction is usually set to 0.9 [9j. The 
dependence of K NLO from the Wilson coefficients C t and C[ at NLO is given by [[]]: 

K nlo (6,z) = Zi< j= 2,7,8 k$>(6, Z ) {Re[Cf(/i b )Cj 0) *^)] +(C W -C5 J .)} + 

k$(6, z) {Re[Cf (^)Cf *(/p,)] + (C 7 - C')} . (16) 

In the previous expression C'j 0) and refer respectively to the LO and NLO contributions 
to the Wilson coefficients CJ t defined as: 


CM) = C!°>M + + 0(a,a 2 ,) 


(17) 


As in the following we are deriving only one-loop formulas for the Wilson coefficients 


C’j'l, Ci = C-°K We will briefly discuss the effects of including Cj l> in section [O. The 


coefficients k i j{5 1 z) used in the our analysis are calculated for 5 = 0.9 and yfz = 0.22 using 
the formulas derived in B0 .The LO branching ratio expression can be easily derived from 
eq. (|T6l) setting k 7 ° 7 = 1 and all the other = 0, giving: 


K L0 = \CM )\ 2 + \c'M)\- 


(18) 


independently of the choice of 5 and z. 


3 CV ; 8 and Cj 8 gluino contributions to b —> sj 

In the following we will focus on the gluino contribution to the Wilson coefficients C’ 7 ,8 and 
Cy 8 . There is only one gluino diagram that contributes to C 7 and C 7 , with the external 
photon line attached to the down-squark line, while two diagrams can contribute to the 
C% and C'g coefficients, as the gluon external line can be attached to the squark or the 
gluino lines. The one-loop gluino contributions to the C 7 , 8 and C' 7 8 coefficients are given 
respectively by: 
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( 20 ) 
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in which x A = rh^/rh 2 ^, with rhg the gluino mass, and m A the mass of the A-th down 
squark eigenstate. L d and R d are the Left and Right gluino couplings to a generic down 
quark d given by: 


L 9 d = -V2U A4 , Rl = V2U A4+3 , (23) 


in which U is the 6 x 6 down-squark rotation matrix. The loop integrals F V2 and F 43 are 
defined as: 


F 21 = F 2 {x) + 9Fi(x) , F 43 = Ft(x) + 9F 3 (x) , 


(24) 


191 for an easier connection with the 


using the conventions for the integrals Fi(x) as in 
standard convention in the literature. 

It is illustrative to write the gluino contribution to the 67,8 and C' 7 8 Wilson coefficients 
using the MI approximation. First, note that the set of integrals used in [|T9| is not the 


most appropriate for dealing with the MI formulas. However, for the sake of simplicity 
we will retain these conventions and further define the integrals iq and their “derivatives” 
through the following self-consistent relations: 


Fi(-) = -U(x,y) , 
y y 


F^{-) = ~fi(x, V ) , 
y y Oy 
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1 1 d n 

n\ y n+1 dy n 
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Using this notation, the hrst and second order terms in the MI expansion for the C 7 g and 
C 7 8 coefficients are given respectively by: 
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In the previous formulas x 9 D = m^/rh 2 Dl with rhr, the average down-squark mass related 
to the down-squark mass eigenstates via the relation m 2 A = fh 2 D + Sm 2 A . The dehnitions of 
the MI parameters are: 


S LL = 

*7 
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E Ut A 6m\U A j 


D A=1 


S RR = 
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E Ul 3 Jm\U AJ+ 3 , 


D A=1 
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(33) 


D ^=1 


In deriving eqs.([2^-^2|) to the second order in the MI parameters, we have kept only the 
dominant term proportional to tan /3 (the A b term is retained in the above expression for 
defining our convention for the /a term; see later), and neglected all of the other off-diagonal 
mass insertions. Clearly the dominant terms in eqs. ( [25| - |32l) are those proportional to the 
gluino chirality flip, such that the gluino contribution to C^ {C' 7 ) at first order depends 
only on the MI term (5 R p ). However, for large tan/1 and /r ss in a, the second order MI 
terms in eqs. ( p9|]32l ) can become comparable in size with the hrst order mass insertions. 
Thus, two different MI parameters are relevant in the L/R sectors: (<5^ 3 R , b R p) and (S R p , 
S R , R ), contrary to common wisdom. To which extent the LL and RR Mis are relevant 
depends of course on the values chosen for /1 and tan (3 , but in a large part of the allowed 
SUSY parameter space they cannot in general be neglected. Moreover, the fact that the 
gluino Wilson coefficients depend on two different MI parameters will have important 
consequences in the study of the b —> sy CP asymmetry^ 


4 Alternative solution to b —>■ S7 branching ratio 


In the majority of the previous studies of the b —> sy process, the main focus was to 
calculate the SM and NP contributions to the g coefficients. The contributions to 
b —> sy coming from C' 7 8 have usually been neglected on the assumption that they are 
suppressed compared to CV ]8 by the ratio m s /m b . While this assumption is always valid 
for the SM and for the Higgs-sector contributions, in the case of the uMSSM this is not 
generally the case. It is only within specific MSSM scenarios (such as MFV) that the 
gluino and chargino contributions to the C 78 coefficients can be neglected due to the 
m s /mb suppression factor. In the general uMSSM this suppression can be absent and, 
in particular, the gluino contributions to CVy and C 7 8 are naturally of the same order 

mm- 

Therefore, in the following we present an alternative approach to the b —► sq process in 
supersymmetric models. We assume a particular scenario in which the total contribution 


6 Specifically, if only the first order term in the MI is taken the b 
discussed in greater detail in section 4.3. 


sy CP asymmetry vanishes, as 
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to C* 7 i 8 is negligible and the main contribution to the b —> sy branching ratio is given by 
Cyg. This “Cy dominated” scenario is realized when the chargino, neutralino, and gluino 
contributions to CV, 8 sum up in such a way as to cancel the W and Higgs contributions 
almost completely^. In our opinion this situation does not require substantially more 
fine tuning than what is required in the usual MFV scenario, where conversely the NP 
contributions to CV, 8 essentially cancel between themselves (or are almost decoupled) so 
that all the measured b —> sy branching ratio is produced by the W diagram. As previously 
discussed, many classes of SUSY breaking models lead to off-diagonal LR entries of the 
down-squark sector that are naturally suppressed compared with those of the up-squark 


sector: 


(O 


max(mj, rrij) 


K R ) U 


m t 


y 


(34) 


in which rriij are down-quark masses. In particular, the entries, which are relevant 

for the b —> sy process, receive a 0 (mb/m t ) suppression as can be derived from eqs.(j^-|j). 
For ~ 0(1), a natural value ) d ~ 0{mb/m t ) ~ ICC 2 is obtained. With this 

mechanism at work, off-diagonal chargino and gluino contributions to flavor changing pro¬ 
cesses are naturally of the same order. The a s /a w enhancement of the gluino contribution 
with respect to the chargino one is compensated by the mb/m t suppression of the LR 
off-diagonal entries. Clearly a complete analysis of the regions of uMSSM parameter space 
where the CV, 8 cancellation takes place is an important task, necessary for studying the 
details of this scenario. However, a detailed analysis is beyond the scope of this paper will 
be discussed elsewhere ]39|]. It is worth mentioning that in preliminary scans we checked 


that it is not difficult to find a candidate set of parameters where CV, 8 numerically yield 
small contributions to the b —> sy branching ratio. Of course, this set is not obviously 
not expected to be unique, and further checking that any such parameter sets are consis¬ 
tent with all the other existing measurements of FCNC and CP-violating observables will 
impose further strong constraints. 

Finally, we stress that in the following analysis we do not make any specific assumptions 
as to the size of the off-diagonal entries of the down-squark mass matrix. In particular, 
we are not using any of the relations described in eqs. (Q-^j). The previous arguments have 
been intended as a theoretical framework for the following model independent analysis. A 
general discussion of the CP-violating sector, using the factorization a n sat z of eq. ([]]), will 
be the subject of a forthcoming paper J39| . 


7 The main constraint on this scenario is the requirements of the C’j cancellation. The Cg, contribution 
enters in the b —> sy branching ratio at 0(a s ) and usually cannot account for more than 10% of the 
measured branching ratio. 
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Figure 1: The dependence ofb^-s'y branching ratio on 8§f and 5^ for different values of 
rrig/nio, for tan (3 = 20 and /i = 350 GeV. All of the other off-diagonal entries except the 
one displayed in the axes, are assumed to vanish. C 7t8 (pw) = 0 is assumed. The horizontal 
lines represent the la experimental allowed region. 


4.1 Single MI dominance analysis 


From eqs. (|25| - P2| ), one can read (in MI language) the off-diagonal entries that are relevant 
for the gluino contribution to the C 7 $ and C' 78 Wilson coefficients^. Note that limits on 
8 f^ ~ G( 10~ 2 ) have previously been obtained in |fl3H - No stringent bound has been derived 
there for 8 ff , as this term at lowest order does not come with the frig/mb enhancement 
(see eqs.(|5|)). No limits were derived on 8 ff 8 and 8 ^ because in the specific scenario used 
in [|13j, the “opposite chirality” Mis are suppressed by a factor m s /m b and so negligible. An 
analysis of the 8 ff 8 dependence has been performed in [^] , in which the W contribution to 
C 7i8 was not set to zero (sometimes also Higgs and MFV chargino contributions to C 7t8 were 
included). Consequently their bounds on the down-squark off-diagonal Mis contributing 
to C 7 8 are more stringent than the bounds we derive in our scenario, for which the total 
contribution to C 7 $ is assumed to be negligible. It is clearly only in the scenario we study 
that an absolute constraint on these Mis be derived. Moreover no analysis on 8 ff and 
8 iff 1 was performed in |28|] as these contributions are not relevant in the small tan f3 region, 
as can be seen from eqs.(^-|32|). 

In Fig. [T] we show the dependence of the b —> sy branching ratio on the MI terms 
82 8 and 8 ^ for different values of x 8 9 D = fh 2 -/fh 2 D and for tan/? = 20 and pi = 350 GeV. 

8 From now on, for the sake of simplicity the symbol Sf B will be used instead of {8f B ) d for referring 

to the down-squark Mis. 
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Figure 2: Dependence ofb^s'y branching ratio on 5Iff* for different three values of tan (3 , 
with the other parameters fixed to fhg/fhq = 350/500 and /i = 350 GeV. All of the other 
off-diagonal entries, except the one displayed in the axes, are assumed to vanish. The 
horizontal lines represent the la experimental allowed region. 

All the other off-diagonal entries in the down-squark mass matrix are assumed to vanish 
for simplicity. “Individual” limits 8 f/f < ICO 2 and S^ < 1.5 x ICO 1 can be obtained 
respectively from the left and right side plot of Fig. [I]. Horizontal full lines represent 
lcr deviations from the experimental results reported in eq. (jTJ) . Of course, the required 
cancellation of the total CYs contribution may in general need nonvanishing off-diagonal 
entries of the up and down squark mass matrices. However, the specific values of these 
entries do not significantly affect the absolute limits on 8 ^ and 8 Iff 1 Mis shown in Fig. |1]. 

As expected from eqs.(|3l|,|32D, the bounds obtained for Htff 1 are strongly dependent 
on the product /Man (3. In Fig. [2| we show the tan/3 dependence of this limit, for fixed 
fhg/fhq = 350/500 and p = 350 GeV. More stringent bounds on 5 /,| R can be obtained 
for larger tan/3 . For tan/3 > 35 the bounds on 8 fff { can become as stringent as the 
Siftf bounds. Similar considerations and bounds obviously hold also for the 8 ft MI. As we 
are only interested here in the gluino contributions to C/ 8 , we do not discuss this sector 
in detail. Clearly this term must be taken into consideration if a similar analysis was 
performed for the CYs coefficient in the large tan /3 region. 

In Fig. [I] and Fig. we set Cj = C$ = 0 so that the only contribution to the 
b —> sy branching ratio is due to the gluino contribution to C’ 7 and C 8 . Thus one should 
think that for vanishing 8 ^ and/or Slff 1 the branching ratio in our scenario should van¬ 
ish. The reason for the finite, nonzero contribution is the fact that we are using a NLO 
formula for the b —> sy branching ratio |J. At NLO, imposing the condition = 0 

still leaves constant terms that arise from the mixing of the SM operators (specifically, 
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Figure 3: Dependence ofb^-s'y branching ratio on 5 ^| L forrhg/rhq = 350/500, tan/? = 20 
and /i = 350 GeV. All the other off-diagonal entries, except the one displayed in the axes, 
are assumed to vanish. In the plots we show the result obtained using LO (dashed line) and 
NLO (full line) formula for the b —► sy branching ratio, setting respectively C 7 ffpw) = 0 
(left plot) and C 7t s(pb) = 0 (right plot). The horizontal lines represent the la experimental 
allowed region. 

in our chosen basis, C-ff that do not contribute to the branching ratio at LO. In Fig. |3] 
(left side), we compare the results obtained using the LO and NLO expression for the 
b —> sy branching ratio imposing the condition C 7t s{^w) — 0. As can be seen explicitly, the 
difference in using the LO or NLO is sizeable. In Fig. [3] (right side), we compare the results 
obtained using the LO and NLO expression for the b —> sy branching ratio imposing the 
condition CV i8 (//*,) = 0. As one can see now, the LO contribution to the b —> sy branching 
ratio vanishes for vanishing Mis. This does not happen for the LO contribution of the left 
plot, as a finite contribution to the branching ratio appears from the running nw —> Hb 
when the condition C 7 ffjiw) = 0 is taken. In all the plots, except Fig. [| (right side), we 
use C 7 ) 8 (/iv y) — 0 ) as this is the natural scale where cancellations could be explained in 
terms of the underlying fundamental theory, while the choice CV >8 (//*,) = 0 seems highly 
accidental. Finally, it should be noted that the strongest restriction comes from imposing 
the condition C 7 = 0. The same requirement on C& could easily be relaxed, and our re¬ 
sults would remain almost unchanged. The C 8 contribution to the b —> sy branching ratio 
represents in fact only a 10% effect of the total. 

It is important to notice at this point that a consistent analysis of b —> sy at NLO would 
require the calculation of the two-loop (QCD and SQCD) contribution to the C 7 coefficient. 
In the general uMSSM the calculation of the O(o%) contribution to C 7 (and obviously C 7 ) 
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is extremely complicated. In |57J , the contribution to C 7 from the two-loop diagrams with 
one gluino and one gluon internal line has been calculated. This represents the dominant 
MSSM two-loop contribution only in the limit of very heavy gluino mass of O(lTeV) and 
small tan/3 (~ 1). Thus it cannot be applied to our analysis, in which SUSY masses (and 
the gluino mass in particular) below 500 GeV and large tan/3 are assumed. In fact, if 
the gluino mass is light the two-loop diagrams with two gluino internal lines should also 
be taken into account. Moreover, if tan/3 is large, diagrams with internal Higgsino lines 
cannot be neglected anymore as Yukawa couplings can become of 0(1). Using the results 
of [[U]], one obtains an effect of a few percent in the b —> sy branching ratio. It should 


be remembered, however, that in our analysis this provides only a very crude estimation. 
It seems reasonable to expect a possible 10% modification of the b —> sy branching ratio 
results from the inclusion of the complete NLO calculation of the C 7 coefficient. Moreover, 
while the two-loop diagrams with gluino/gluon internal lines have the same MI structure 
and as such are proportional to the one-loop gluino contribution to C 7 , this is not the case 
for the diagrams with gluino/Higgsino internal lines, for which the CKM flavor changing 
structure also enters. 


4.2 General MI analysis 

A general analysis of the gluino contribution to C 78 depends simultaneously on both 
the 5 23 and iR Mis. For a complete specification of our scenario the only other free 
parameters that need to be fixed are the ratio between the gluino mass and the common 
down-squark mass, 7hg/m D , the product /j, tan/3 , and the relative phase between 5 and 
8 RR . The influence in of all the other down-sector squark matrix off-diagonal entries 
and MSSM parameters in the C 7a sector can safely be neglected^. Thus, we can have a 
complete description in terms of only five free parameters of the b —> sy phenomenology 
in our MSSM U C 7 dominated” scenario. 

In Fig. f| we show the la experimentally allowed region in the (8 R 8 , 8 ^ R ) parameter 
space for a specific choice of frig/friD = 350/500, /i = 350 GeV, and for three different 
values of tan/3 =3, 20 and 35. For 5 RR or S RR vanishing, one obtains the regions depicted 
in Figs. HI and |2j. Larger regions in the (8%$ , 8 RR ) parameter space are obtained when 
both the Mis take nonvanishing values. It is clear no absolute limit can be derived for 
the two Mis simultaneously. The values (8 RR , 8 RR ) ~ (1,0.1) are, for example, possible^ 
for tan/3 = 35. In fact, as can be seen in Fig. [|, there is always a “flat direction” where 

9 Of course all of the other off-diagonal entries of the down-squark and up-squark mass matrices as well 
all the other flavor conserving MSSM parameters enter in our analysis, as we assume to choose them in 
such a way that the condition C 7} s = 0 is satisfied. However, as previously mentioned the detailed analysis 
of this condition will be discussed in a following paper Q. 

10 One should check if, for such large MI values, charge and color breaking minima appear. Anyway as 
these are usually rather model depend assumptions we don’t introduce here the constraints discussed, for 
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Figure 4: la-allowed region in the (Sff[ { , <5|| L ) parameter space for three different values 
oftan/3, with the other parameters fixed to fhg/fhq = 350/500, and /i = 350 GeV. All the 
other off-diagonal entries, except the one displayed in the axes, are assumed to vanish. 


large values of 8 ffp and Sffp can be tuned in such a way that the gluino contribution to 
C 7 8 is consistent with the experimental bound. This flat direction clearly depends on the 
chosen values for fhg/fhq and /it an fj . The presence of this particular direction is explained 
by the fact that we are allowing complex off-diagonal entries. Hence the relative phase 
between 8 §f and Sf/p can be fixed in such a way that the needed amount of cancellation 
can be obtained between the first and second order MI contribution. In the notation used 
in eqs. (|27| , |3T|) the line of maximal cancellation is obtained for <p = arg[h|| i 5§P } = ±7T. 


example, in ft^. 
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4.3 CP asymmetry and branching ratio 

In addition to the b —> sy branching ratio, the experimental collaborations will provide in 
the following years more precise measurements of the b —» 57 CP asymmetry: 


A C p(b sy ) 


BR{b —> 57 ) — —> sy ) 

> sy ) + —> 57 ) 


(35) 


The present best experimental value available [41 ] gives at 90% CL level the following 
range: 

- 0.27 < A C p{b -> sy ) < 0.10 , (36) 


which is still too imprecise for to provide useful tests for NP, although the measurement 
is expected to be upgraded soon. 

The only flavor-violating and CP-violating source in the SM (and MFV scenarios) is 
given by the CKM matrix, which results in a very small prediction for the CP asymmetry. 
In the SM an asymmetry approximatively of 0.5% is expected [9;j. If other sources of CP 
violation are present, a much bigger CP asymmetry could be produced (see references || 
and ||42| ). 

I 11 our C 7 dominated scenario, one can derive the following approximate relation for 
the CP asymmetry [9j], in terms of the | L and S.™ Mis: 


ah \ 4 •ImtCyC 

A C p(b -> sy) = --a s (fi b )— 


k(x 9 D 


f iribii tan /3 

V ™ 2 d 


Isin p , (37) 


in which ip is the relative phase between and <5|| R as previously defined. The constant 

of proportionality k(x 9 D ) depends only on the ratio fhg/fhr, through the integrals F t and can 
be easily obtained from eqs.(|2^j3l|). One can immediately note from eq.(]37|) that if only 


one MI is considered, the CP asymmetry is automatically zero. Moreover, a nonvanishing 
phase in the off-diagonal down-squark mass matrix is necessary^]. No sensitive bounds on 
this phase can be extracted from EDM’s in a general flavor violating scenario. 

I 11 Fig. |j. we show the results obtained for the branching ratio and CP asymmetry in 
which 5^ , S™ and the relative phase <p are varied arbitrarily for a fixed value frig/frig = 
350/500 and tan/3 = 35. The full vertical lines represents the lcr region experimentally 
allowed by the b —> sy branching ratio measurements. It is possible, using C' 7 8 alone, to 
saturate the b —*■ sy measured branching ratio and at the same time have a CP asymmetry 
even larger than ± 10 %, the sign of the asymmetry being determined by the sign of sin<p. 
As Fig. H shows, in the relevant branching ratio range the CP asymmetry range is constant. 
No strong dependence from tan/3, in the large tan/3 region, is present. The points with 


11 Recall that for avoiding EDM constraints reparameterization invariant combinations of flavor- 
independent phases (such as the phase of p in a particular basis) are taken to be zero. 
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Figure 5: Asymmetry vs branching ratio for three different values of tan fj, with rhg/rhq = 
350/500, and /i = 350 GeV. All the off-diagonal entries except | L and 5™ are assumed 
to vanish. The vertical lines represent the la experimental allowed region. 


large asymmetry (> 5%) lie in the “flat direction” observed in Fig|4| and they have almost 
p ps ±7 r (obviously for p = ±7r the CP asymmetry vanishes). The explanation of this fact 
is the following. The numerator is proportional to sin<^ and so goes to 0 as <p approaches 
± 7 T. However, at the same time it is enhanced for large MI values. This happens when 
the flat direction condition is (almost) satisfied. Here, in fact, a cancellation between the 
two (large) MI terms takes place, providing the enhancement of the CP asymmetry as the 
denominator remains practically constant, fixed by the allowed experimental measurement 
on the branching ratio. Note also that for parameter values outside the flat direction 
condition a CP asymmetry of a few % can still be observed, about ten times bigger than 
the SM prediction. The same order of magnitude can be observed in MFV when large 
tan f3 effects are taken into account . In our scenario even smaller values of the CP 
asymmetry can be obtained, e.g. if one of the two off-diagonal entries is negligible, or the 
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two Mis are “aligned”. 


4.4 Distinguishing the U C' 7 dominated” scenario from the “C 7 
dominated” 


A possible method for disentangling the relative contributions to the b —> s 7 branching 
ratio from the Q 7 and Q\ operators utilizes an analysis of the photon polarization. A 
detailed analysis of how it is possible to extract information from the photon polarization 
in radiative B decays is given in JJ3J. For simplicity, let us define the following “theoretical” 
LR asymmetry at LO: 


A L ii{b -*■ 57 ) = 


BR{b - s lL ) - BR(b - s^ R ) \C 7 (^)\ 2 - \C'M? 


BR(b -► S7l ) + BR{b -► s lR ) |C 7 (^)| 2 + | C 7 {ii b )\- 


(38) 


which could in principle disinguish between C 7 or C 7 dominated scenarios. Here L,R is 
the polarization of the external photon. This quantity is related to the quark chiralities 
of the Q 7 , Q' 7 operators. Note that the photon polarization is the best possibility to gain 
information on the operator chirality, which gets almost lost in b and s quark hadronization 
into spin zero mesons (in principle if hadronization into spin one states could be isolated, 
perhaps some information could be obtained). Such a measurement is not yet available as 
only the average quantity BR{b —> s 7i ) + BR{b —> s^ R ) is reported experimentally. 

I 11 the SM case, and in general in all the MFV and mSUGRA scenarios, only the 
C 7 coefficient gives a nonnegligible contribution to the b —> s 7 branching ratio. Only the 
right-handed bottom quark (in the center of mass reference frame) can decay, producing 
a photon with Left polarization and Auiib —> s 7 ) = 1. Small deviations from unity are 
possible due to subleading m s /mb terms and hadronization effects. In our scenario, where 
the C 7 contribution is negligible, only left-handed bottom quarks can decay, emitting a 
photon with Right polarization, which in turn predict Alr{1> —> s 7 ) = —1. In any other 
MSSM scenario, with nonminimal flavor violation, any LR asymmetry between 1 and — 1 
is allowed. Consequently, a measurement of Ai R (b —> s 7 ) different from 1 will be a clear 
indication of physics beyond the SM with a nonminimal flavor structure. It will be very 
interesting to know if (and how precisely) CLEO, BABAR, and BELLE can measure the 
LR asymmetry of eq. (p 8 |) . 


5 Conclusions 

In this letter, we have discussed an alternative explanation of the b —> s 7 branching ratio 
in the MSSM with a nonminimal flavor structure. We analyzed in particular the gluino 
contribution to the Wilson coefficient C 7 associated with the “wrong” chirality operator 


19 





Q' 7 . We show that this coefficient arises mainly from two off-diagonal entries: 5 ^| L and ■ 
For scenarios in which where the C 7 $ contributions to b —> sy are small, (i.e. for regions 
in the MSSM parameter space where W, Higgs, chargino and gluino contributions to CVy 
tend to cancel each other) C 7 8 provides the dominant effect. We derived absolute bounds 
separately on each of these coefficients. We then described the allowed region of (S ^, 
S.™) parameter space, as a function of tan/?. We observed that (for a fixed ratio rhg/fhq 
and for each chosen value of n tan/?), there exists a “flat direction” where large (even 
0(1)) off-diagonal entries are allowed. Along this direction the relative phase between the 
two MI elements is tp ~ tt. For the majority of parameter space in this scenario the CP 
asymmetry is less than 5%. Asymmetries as big as 20% can be obtained along the “flat 
directions”. Finally, we suggested a possible quantity (a LR asymmetry) that (if measured) 
can help to disentangle the C 7 from the C 7 contribution to the b —> sy branching ratio. Any 
Alr(J> — 5 ► sy ) % 1 would be an irrefutable proof of physics beyond the SM. In addition, 
in the framework of the general MSSM, it would indicate the existence of nonminimal 
flavor violation produced by off-diagonal entries in the down-squark mass matrix, generally 
related to a nonzero gluino contribution. In our U C 7 dominated” scenario, where the gluino 
contribution produce the only “visible” effect, we obtain in particular the extreme value 
A.Ln(b —> sy ) = — 1. It would be very interesting if such a quantity could be measured. 
One implication of our analysis is that previous results on MSSM parameters, including 
constraints on the “sign of /i” (i.e. its phase relative to A t ), are more model dependent 
than have been generally assumed. 
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